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Abstract 

We present the resuhs of a complete tree level calculation of the 
processes pp{p) — > Wbb and Wbb + jet that includes the single top 
signal and all irreducible backgrounds simultaneously. In order to 
probe the structure of the Wtb coupling with the highest possible 
accuracy and to look for possible deviations from standard model pre- 
dictions, we identify sensitive observables and propose an optimal set 
of cuts which minimizes the background compared to the signal. At 
the LHC, the single top and the single anti-top rates are different 
and the corresponding asymmetry yields additional information. The 
analysis shows that the sensitivity for anomalous couplings will be im- 
proved at the LHC by a factor of 2-3 compared to the expectations for 
the first measurements at the upgraded Tevatron. Still, the bounds 
on anomalous couplings obtained at hadron colliders will remain 2-8 
times larger than those from high energy 76 colliders, which will, how- 
ever, not be available for some time. All basic calculations have been 
carried out using the computer package CompHEP. The known NLO 
corrections to the single top rate have been taken into account. 



1 Introduction 



The observation by the CDF and DO collaborations of a very heavy top 
quark with a mass of about 175 GeV, close to the indirect prediction from 
fits of precision electroweak data 177t?t}^ GeV 0, has been an important 
confirmation of the Standard Model (SM). Still, important open problems 
remain: why is the top is so heavy and is it really a point like particle? 
A curious numerical coincidence between the mass and vacuum expectation 
value v/\/2 = 175 GeV sets the top quark Yukawa coupling close to unity. 
As has been stressed some time ago 0, because of such unique properties, 
the top quark might provide for the first time a window to the physics of 
electroweak symmetry breaking and even "new physics" beyond. 

A possible signal would be deviations from the SM predictions of the 
interactions of the top quark with other fields. Therefore it is important 
to study and measure all top quark couplings, in particular the coupling to 
the VT-boson and the 6-quark which is responsible for almost all top quark 
decays. Therefore, events with the production of a single top quark are 
extremely interesting at different colliders, because they are directly pro- 
portional to the Wtb vertex. Thus one can hope to measure the structure 
of the vertex and possible deviations from the SM predictions with a high 
accuracy. Furthermore, it should be noted that the processes of single top 
production and decay involve light Higgs production simultaneously (cf. the 
discussion in [Q). The measurement of the Wtb vertex in 76 collisions has 
been described in 0, ||. 

In this paper we discuss the possible accuracy in the determination of 
the structure of the Wtb vertex at the upgraded Tevatron collider and at the 
LHC. Improving on previous considerations we perform a complete tree 
level calculation, taking into account contributions from anomalous operators 
to the Wbt vertex, the production of Wbb and Wbb+jet, which includes the 
single top signal together with the irreducible backgrounds. We have included 
the NLO corrections to the single top part [0]. Based on an analysis of the 
singularities of Feynman diagrams and on explicit calculations we identify 
the set of the most sensitive variables and their corresponding optimal cuts. 
This allows us to obtain a clean single top sample above the background and 
a handle on possible deviations from the SM expectations. 
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2 The Basic Processes 



Single top production at hadron colliders has been studied by a number of 
authors (cf. H, 0, |H and references therein). So far, the most complete set of 
SM processes contributing to the single top rate has been studied in p and 
the most accurate NLO calculations to the main processes have been pre- 
sented in [0]. In recent papers 0, complete Monte Carlo (MC) analysises 
of the single top signal versus backgrounds have been presented. The Feyn- 
man diagrams for all processes contributing to the single top production rate 
have been presented previously (cf. e. g. 0) and include virtual W s-channel 
exchange, ly-gluon fusion and W + top production. The first process is the 
simplest 2^2 reaction, while the VT-gluon process includes 2 — >■ 3 parton 
diagrams. In order to resum the large QCD corrections from g ^ bb splitting 
in the latter, they are combined with the 2 — 2 process involving a b quark 
in the initial state |^, |^ and the corresponding piece of the g ^ bb splitting 
function is subtracted in order to avoid double counting. 

Finally, the W + top production process gives a large contribution to 
total single top production at the LHC ^ . However we will not consider this 
process here, because it does not contribute to the topologies which we will be 
interested in. Moreover, it has a final state similar to top pair production and 
after a suitable background subtraction it is less sensitive to the Wtb vertex 



(cf. the discussion of the corresponding process in e~^e~ collisions |TT[). 

Irrespective of the different strategies employed, the analysises of 
and fl^ have both shown that the single top production rate is large enough 



to be visible above the backgrounds using proper cuts. This is possible despite 
the fact that the background reduction is much more complicated compared 
to the case of top pair production. 

However, in order to probe the Wtb vertex, one must check for possible 
deviations from the SM predictions and then all SM contributions become 
part of the background. Therefore it is necessary to find even stronger cuts, 
defining phase space regions where the deviations from the SM predictions 
for single top production will be most prominent. Approaching this problem 
we perform an accurate calculation of the two processes 

pp — > bbW (la) 

and 

pp^bbW+jet (lb) 
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Figure 1: Feynman diagrams for the process ud hhW . 

which include simultaneously the single top signal and the irreducible back- 
grounds. The Feynman diagrams for process ( p^aj ) at parton level are shown 
in figure |I| and the Feynman diagrams for process ud — *• hbW + jet as a rep- 
resentative of the processes contributing to the process ([T^) at parton level 
are shown in figure |^. Parton processes with gluons in the initial state are 
not shown, but are included in the calculation. The diagrams include the top 
signal, Higgs contribution, QCD diagrams with gluons in the intermediate 
state and several other electroweak diagrams which have to be taken into 
account for an accurate calculation of the rate of very hard processes. Here, 
the Higgs contribution is part of the background to the single top production. 
In our calculation, we have assumed a light Higgs with a mass in the range 
80-120 GeV as a worst case scenario for the single top signal, since the Higgs 
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Figure 2: Feynman diagrams for the process ud — > ghhW . 
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rate drops rapidly with rising Higgs mass \T7]. Even in this case, the Higgs 
contributions will turn out to very small in the phase space regions which 
will be important for single top production. 

All the calculations have been performed with the computer program 
CompHEP , including the proper mapping of singularities and smoothing 
in singular variables in 



We have used the NLO CTEQ4M parametrization 
of parton distribution functions [|15]. For the quark induced processes (|IaD, 
we have chosen the QCD factorization scale to be Mt. This choice is dic- 
tated by the fact that we are selecting a kinematical region where the two 
quarks annihilate into a state close to the top quark mass shell. For pro- 
cesses involving VT-gluon fusion ([TEI) , the choice of scale is more subtle, as 
has been pointed out by [|^ . Therefore, we have pragmatically fixed the scale 
by matching our LO cross sections to the NLO results of 0. This procedure 
leads us to a factorization scale of ~ (Mi/2)^. The fact that this scale is 
reasonably close to the top quark mass shows that the corrections are not very 
large and serves as a a posteriori justification of the pragmatical procedure. 
Therefore we have taken into account the important parts of the NLO cor- 
rections in the hard kinematical region which we are interested in. Finally, 
we notice that the requirement of a double 6-tag in the hard kinematical 
region under consideration suppresses the contributions from the processes 
with the 6-quark in the initial state. Therefore this source of theoretical 
uncertainties for the signal is absent in this case. 



3 Anomalous Wtb Couplings 



In the model independent effective Lagrangian approach [jT6l seven anomalous 
CP conserving operators of dimension six contribute to the Wtb vertex with 
four independent form factors (cf. for explicit expressions). In this paper, 
we do not attempt a simultaneous analysis of all seven operators. Instead we 
study two anomalous operators of the magnetic type as an example in order 
to explore the potential of the colliders. In fact, the V — A coupling is as 
in the SM with the coupling Vtb very close to unity, as required by present 
data . A possible V + A form factor is severely constrained by the 
CLEO b s'y data on a level which is stronger than expected even at 
high energy 76 colliders 0] . This leaves us with the remaining two magnetic 
form factors and we have studied them in the processes 

As before Ef, we have adopted the notation for the Lagrangian in unitary 
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gauge from []TH 
v/2 



2M, 



w 



+ h.c. 



(2) 



with iy± = D^W^ - D,W^, D^ = d^- leA^, a^^ = z/2[j^, 7,] and P± = 
(1 ± 75)/2. The couphngs F2 and are proportional to the coefficients 
Ctw<i> and C{,iy$ of the effective Lagrangian 



.L(i?) _ Ct{b)w^ \/2vM] 



w 



A2 



(3) 



The resulting Feynman rules (cf. the appendix of (Q) have been implemented 
in CompHEP. 



4 Sensitive Variables, Background Suppres- 
sion and the Structure of Singularities in 
Feynman Diagrams 

As always, the correct mapping of the singularities of the Feynman diagrams 
is absolutely crucial in order to achieve numerically stable results in the MC 
integration over phase space. In this section, however, we will focus on a 
different aspect of the singularities, that has impact on physics analysises. 

Unless particular cuts are applied, most of the contributions to the rate 
of any given process come from phase space regions close to the singularities. 
Indeed, this simple observation forms the basis for most of our intuition for 
finding selection cuts to enhance the signal. In many cases, it is however also 
possible to reverse this argument and use the singularities to systematically 
devise cuts for background suppression. This approach requires an analysis 
of all Feynman diagrams contributing to the process, signal and background. 

Shifting the focus from signal enhancement to background suppression 
in this manner is useful for processes with a high rate of both signal and 
background that can afford to lose some rate. We will demonstrate that 
single top production falls into this class. 

The general procedure compares the set S of variables with singularities 
from all signal diagrams with the same set B from all background diagrams, 
reducible and irreducible. If B\S ^, i.e. there are variables with singulari- 
ties in the background diagrams which are regular in all signal diagrams, then 
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it is obvious that singularities in these variables should be cut out as strongly 
as possible. It turns out that the number of different singular variables is 
very limited in the cases of practical interest and a general classification al- 
lows recommendations for choosing sensitive variables. The application of 
this approach to neural net methods will be discussed in . 



Shifting our attention back to the special case of single top production, we 
note that the single signal diagram for ([l^) has only one singular variable: the 
invariant mass Mwb of the top decay products near the top pole Mwb = Mf. 
These contributions have to be kept, of course. 

In the background diagrams, the only s-channel singularities are in the in- 
variant mass Mfji of the bb pair at 0, Mz, and Mh from the coupling to neutral 
vector bosons and Higgs particles. Since the CKM matrix element Vub is tiny, 
we can ignore the multiperipheral diagrams with the Wub-conpling and we 
only have to consider the t-channel variables t^^fefe = tj-^vK td^bh = tu^w- 
Unfortunately, the t-channel variables are not directly observable in hadron 
collisions. We can, however, use the corresponding transverse momentum as 
a surrogate. In this case, these are the P* of the 66-pair or, equivalently, of 
the M^-boson. 

From this simple argument, we conclude that the invariant mass M^j and 
the transverse momentum P^, are the most effective variables for expressing 
cuts for the process (|T^. 

Analogous considerations for the diagrams in figure |^ leads to the vari- 
ables M^i, P^j^ and (the latter are now no longer equivalent) for the 
process (pIB]). Here, the transverse momentum distributions of single jets Pi 
and P* are problematic, because the same singularity occurs for both signal 
and background diagrams and the signal and the background will have simi- 
lar shapes. Therefore cuts on these variables must be defined with discretion 
to achieve a balance between the competing goals of good jet identification 
and high signal rate. 

Of course, there are more variables that can have different distributions 
for the signal and background, but which are not directly related to the 
singularities of Feynman diagrams. One such variable is the partonic center 
of mass energy s. The difference is caused here by different thresholds for 
the signal and the backgrounds. 
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Figure 3: Distributions for invariant masses and transverse momenta for the 
process pp — > bbW at Tevatron using the soft cuts in (^). 
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Single Top (x5) and complete Wbb process (Vs = 14 TeV) 
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Figure 4: Distributions for invariant masses and transverse momenta for the 
process pp bbW at LHC using the soft cuts in (^). 
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Single Top(x3) and complete Wbbj process (Vs = 2000 GeV) 
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Figure 5: Distributions for invariant masses and transverse momenta for the 
process pp — » jbbW at Tevatron using the soft cuts in (HH). 
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Figure 6: Distributions for invariant masses and transverse momenta for the 
process pp — » jbbW at LHC using the soft cuts in (^). 
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5 Distributions and numerical results 



To illustrate the kinematical properties of the processes (|l|) we show in fig- 
ures H, H, 1^, and P several distributions on the variables discussed above with 
soft initial cuts on the jet Pj, jet rapidity and jet cone size ^Rjj^ej) 

P/>10GeV, \r]j\<2.h, APjj(ej) > 0.5 } Tevatron (4a) 
P/>20GeV, |r?j|<3, APjj(ej) > 0.5 } LHC (4b) 

The figures allow to compare the distributions from the single top part only 
with those from the complete set of the SM diagrams. The notations hi 
and 62 refer to the 6-jets with the larger and smaller P* respectively. To 
make the contribution of the single top more visible in the figures, we have 
scaled the rate by appropriate factors as indicated. As we have expected from 
the analysis of the singularities, the distributions are significantly different 
for single top and background contributions. 

The additional power of momentum in the amomalous couplings (^ will 
cause a deviation from the SM prediction that rises with energy and transver- 
sal momentum. However, since the rate falls off quickly with P*, the optimal 
cuts must not be too strong in order to conserve rate. 

The optimized cuts turn out to be different for the Tevatron and the LHC, 
as well as for processes with two 6-jets and processes with two 6-jets and one 
light quark or gluon jet on the other hand. For the process ( p!a| ) we find 



Pi > 30 GeV, P4 >20 GeV, 

Mfe5 >100 GeV, P^*5, P^ >30 GeV 

Pb\ > 50 GeV, P4 > 20 GeV, 

M^l >100 GeV, P*^, P^ >100 GeV 



Tevatron (5a) 
LHC (5b) 



and for process ([T^) 



PL >40 GeV, P4 , P* >20 GeV 
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Tevatron (5c) 



M^l >40 GeV, P^ >30 GeV, P^ >20 GeV 

P^ > 50 GeV, Pl,,P-> 20 GeV, 
Mfes > 100 GeV, P^*5 > 100 GeV, P^ >30 GeV 

As we will see below, it is of crucial importance to use both processes (|l£ 
and (llE) for establishing limits on anomalous couplings, in particular at the 



LHC (5d) 
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Process 


Tevatron 
cr/pb 


LHC 
cr/pb 


ud W+hh 
/ ud-^ W-bb 


soft 
cuts (4a) 


optimized 
cuts (|5aD 


sof 
cuts ( 


t 

4b 


) 


optimized 
cuts (^) 


complete 
single top 


8.1 
0.57 


0.68 
0.30 


16.6 / 10.4 
3.2 / 1.8 


3.8 / 2.4 
1.7 / 0.9 


ug dW+bb 
/ ug-^ dW-bb 


soft 
cuts pa|) 


optimized 
cuts (|5c|) 


soft 
cuts 


optimized 
cuts (|^) 


complete 
single top 


1.4 
0.42 


0.32 
0.27 


28.4 / 5.8 
18.0 / 2.0 


9.6 / 1.8 
7.8 / 1.5 


ud gW+bb 
/ ud^ gW~bb 


soft 
cuts (H) 


optimized 
cuts (15^) 


sof 
cuts ( 


t 

4b 


) 


optimized 
cuts (H) 


complete 
single top 


2.5 
0.38 


0.34 
0.13 


4.6 / 1.4 
1.4 / 0.7 


2.6 / 0.8 
0.8 / 0.4 


gd uW+bb 
/ gd^ uW-bb 


soft 
cuts 


optimized 
cuts ([5^) 


soft 
cuts 


optimized 
cuts (0) 


complete 
single top 


0.41 
0.12 


0.08 
0.07 


6.0 / 15.2 
4.0 / 9.0 


1.7 / 4.0 
1.6 / 3.6 



Table 1: Single top cross sections at Tevatron and LHC. The numbers of the 
cuts refer to the formulae in the text. 



LHC. In order to demonstrate the effect of the cuts, we have given the cross 
sections for several sub-processes at both colliders in table |l[ We stress that 
the rates of single top and single anti-top production differ together with their 
corresponding backgrounds at the pp-coUider LHC, while they are equal at 
the pp-collider Tevatron. 

The numbers labeled 'complete' correspond to the contribution from 
all SM diagrams including the single top single and all interferences. The 
numbers show that the optimized cuts @ do indeed improve the signal to 
background ratio significantly. In particular the gluon initiated subprocesses 
provide a clean sample that is dominated by single top production in W^-gluon 
fusion. 

The dependence of the total cross section for the processes (|l|) on anoma- 
lous couplings after optimized cuts (||) is shown in the upper part of figure ^ 
for the Tevatron and in figure |^ for the LHC. The resulting two standard 
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Complete process with anomalous couplings, Vs =2000 GeV 
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Figure 7: Cross sections after optimized cuts (^,^ and corresponding limits 
on anomalous couplings at the Tevatron. 
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Figure 8: Cross sections after optimized cuts (|5b| , |5d| ) and corresponding limits 
on anomalous couplings at the LHC. 
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Systematics 






± 10% 


-0.18 ... 


+0.55 


-0.24 . . . 


+0.25 


± 0% 


-0.07 ... 


+0.11 


-0.18 ... 


+0.21 



Table 2: Uncorrelated limits on anomalous couplings from Tevatron mea- 
surements with and without systematical uncertainties. 



deviation exclusion contours are presented in the lower part of these fig- 
ures. These exclusion contours correspond to the electronic and muonic decay 
modes of the l^-boson, including r cascade decays. The combined selection 
efficiency in the hard kinematical region under consideration, including the 
double 6-tagging, is assumed to be 50% and as integrated luminosities we 
have used 2fb~^ for the upgraded Tevatron and lOOfb"^ for the LHC 

The combined annulus in figure |^ corresponds to the optimistic scenario 
when only statistical errors are taken into account. A systematic uncertainty 
of about 10% is expected for the upgraded Tevatron (cf. the last paper in P). 
The resulting exclusion contour is shown in figure |^ as well (the allowed region 
now covers the hole of the annulus). 

Figure ^ demonstrates that it will be essential to measure both pro- 
cesses (|l]) at the LHC. The allowed regions for the each process alone are 
rather large annuli, but the overlapping region is much smaller and allows 
an improvement of the sensitivity on anomalous couplings by an order of 
magnitude from the Tevatron to the LHC. 

The rate of single top production at LHC is different from the rate of sin- 
gle anti-top production. This asymmetry provides an additional observable 
at LHC that is not available at the Tevatron. The dependence of the asym- 
metry after optimized cuts (pE,|SH|) on anomalous couplings and the resulting 



two standard deviation exclusion contours are shown in figure ^. While the 
allowed region for the process (|T^ is different from the region derived from 
the rate, the combined limits from both processes (|I]) are similar. 

Systematical uncertainties (from AMw, ^Mt, parton distribution func- 
tions, QCD scales, out of cone corrections, luminosity determination, etc.) 
will play an important role at the LHC as well. It is however impossible to 
predict them accurately before cross checks from independent measurements 
at the LHC scale can be performed. Therefore we simply take a set of com- 
bined systematic uncertainty and include it into a new fit for each value. 
Figure ^ shows how the exclusion contours deteriorate when systematic er- 
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Asymmetry with anomalous couplings, Vs =14 TeV 
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Figure 9: Top-anti-top asymmetry after optimized cuts (^,|^) and corre- 
sponding limits on anomalous couplings at the LHC. (Note that the axes are 
exchanged with respect to figures and ||.) 
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Figure 10: Dependence of the limits on anomalous couplings from LHC mea- 
surements on the size of systematical uncertainties. 
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Systematics 


^2 




± 10% 


-0.094 . . . 


+0.34 


-0.17 ... 


+0.18 


± 5% 


-0.052 . . . 


+0.097 


-0.12 ... 


+0.13 


± 1% 


-0.013 ... 


+0.014 


-0.05 ... 


+0.06 


± 0% 


-0.003 . . . 


+0.003 


-0.022 ... 


+0.03 



Table 3: Uncorrelated limits on anomalous couplings from LHC measure- 
ments with for several estimates of systematical uncertainties. 





^2 


Fi< 


Tevatron (A^ys. ^ 10%) 


-0.18 


... +0.55 


-0.24 


... +0.25 


LHC (Asys. ~ 5%) 


-0.052 


. . . +0.097 


-0.12 


... +0.13 


76 (^Se+e- = O.STeV) 


-0.1 


... +0.1 


-0.1 


... +0.1 


76 (Vse+e- = 2.0TeV) 


-0.008 


. . . +0.035 


-0.016 


... +0.016 



Table 4: Uncorrelated limits on anomalous couplings from measurements at 
different machines. 

rors of 1% and 5% are included. In the tables H and ^ for Tevatron and 
LHC respectively, the uncorrelated bounds on the anomalous coupling pa- 
rameters F2 and F.^ are given, assuming different systematic uncertainties. 
Unfortunately, including a 10% systematic error at the LHC will diminish 
the sensitivity significantly and the allowed regions will be comparable to 
those obtained at Tevatron. 

The potential of the hadron colliders should be compared to the potential 
a next generation e~^e~ linear collider (LC) where the best sensitivity could 
be obtained in high energy 7e-collisions [^, |^ . The results of this comparison 
are shown in the Table ^. One can see that a 500GeV LC will outperform the 
Tevatron (assuming a systematic uncertainty of 10%) by a factor of two to 
five. Nevertheless, the upgraded Tevatron is expected start with physics runs 
long before a LC. The upgraded Tevatron will therefore be able to perform 
the first direct measurements of the structure of the Wtb coupling. 

The LHC will only be able to rival a 500GeV LC, when the systematic 
uncertainties can be kept very small (on the order of 1%). This goal will 
be very difficult to achieve. In the more realistic scenario of 5% systematic 
uncertainties, the LHC will improve the Tevatron limits considerably, but it 
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will fall short of a high energy LC by a factor of three to eight, depending 
on the coupling under consideration. 

In the present analysis, we have not included sources of reducible back- 
ground |1^] to single top production at hadron colliders. However, this 
reducible background is sufficiently suppressed in the kinematical regions 
corresponding to our optimized cuts. More detailed simulations and the 
actual analysis should nevertheless include the tails of these background dis- 
tributions as well. 

Finally, we mention that the exclusion contours in figures ^ and |10 
be combined with constraints [0, |] on the right-handed coupling 

< from the CLEO measurement of 6 — > 57 [^] to improve the limit on 
the left-handed coupling F2 ■ 



can can 
0.0015 < 



6 Conclusions 

We have presented the results of a complete tree level calculation of the pro- 
cesses pp{j>) Wbb and Wbb + jet, taking into account the contributions 
of anomalous Wtb operators. These final states simultaneously include the 
single top signal with subsequent decays and irreducible Standard Model 
backgrounds. We have determined the most sensitive variables from an anal- 
ysis of the singularities of the Feynman diagrams in phase space, in order to 
achieve optimal background suppression without sacrificing too much of the 
signal. 

It was shown that the optimized cuts allow to suppress the background 
rate drastically and to extract limits on anomalous coupling parameters. The 
accuracy at the LHC is expected to be better by a factor of two to three com- 
pared to the upgraded Tevatron. Nevertheless, the Tevatron measurements 
will provide the first direct information on the structure of the Wtb coupling. 
For the higher accuracy at LHC, it is essential to measure the two final states 
separately, to perform efficient double 6-tagging at high P* and to control the 
systematic uncertainties at a leverl betetr than 10%. 

At the LHC, one can reduce the dependence of the results on parton dis- 
tribution functions, QCD scales, etc. by using the asymmetry of single top 
and single anti-top production, remains significant after all cuts. Reducible 
backgrounds are expected to be less important in the phase space region 



corresponding to the optimized cuts R ITO[. Nevertheless, a complete sim 



ulation including reducible backgrounds and realistic detector response will 
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be required for the final experimental analysis. 
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